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基于 降 维 Gabor 特征 和 决策 融合 的 高 光谱 图 像 分 类 


杨 秀 杰 !， 高 NP? 
(1. 重庆 电子 工程 职业 学 院 ,数字 媒体 学 院 , 重庆 4013312. 西南 大 学 学 生 工作 处 , 重庆 400715) 


摘 要 : 针对 传统 高 光谱 图 像 分 类 算法 忽略 空间 特征 这 个 问题 ， 提 出 一 种 基于 Gabor 特征 和 决策 融合 的 高 光谱 图 像 
分 类 算法 。 首 先 ， 通 过 系数 相关 甜 阵 智能 地 对 相 邻 和 高 相关 光谱 带 进行 分 组 ; 接着 ， 在 PCA 投影 子 空 间 中 提取 每 组 
中 的 Gabor 特征 ， 以 量化 局 部 方向 和 尺度 特征 ; 然后 ， 结 合 保留 非 负 和 珑 阵 分 解 的 局 部 性 以 减少 这 些 特征 子 空 间 的 维 
度 ; 最 后 ， 对 降 维 特征 进行 高 斯 混合 模型 分 类 ， 并 使 用 对 数 分 类 池 决 策 融 合 规则 将 分 类 结果 合并 。 实 验 结果 表明 ， 
所 提 算 法 优 于 传统 和 现 有 的 共计 8 种 先进 的 分 类 算法 。 
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Hyperspectral image classification based on dimensionality reduction 
Gabor feature and decision Fusion 


Yang Xiujie!, Gao Li? 
(1. School OF Digital Media, Chongqing College of Electronic Engineering, Chongqing 401331, China; 2. Dept. of Student 
Work, Southwest University, Chongqing 400715, China) 


Abstract: Aiming at the problem of ignoring spatial features in traditional hyperspectral image classification algorithm, this 
paper proposed a hyperspectral image classification algorithm based on dimensionality reduction Gabor feature and decision 
fusion. Firstly, the adjacent and hypercorrelated spectral bands are intelligently grouped by coefficient correlation matrix. 
Then, Gabor features in each group are extracted from the PCA projection subspace to quantify the local direction and scale 
features. Then, the dimensionality of these feature subspaces is reduced by preserving the locality of the decomposition of 
nonnegative matrices. Finally, the reduced dimension features are classified by Gaussian mixture model, and the 
classification results are merged by decision fusion rules. Experimental results show that the proposed algorithm is superior 
to eight kinds of traditional and existing advanced classification algorithms. 


Key words: hyperspectral image; classication; Gabor features; Gauss mixed model; decision Fusion; PCA projection 
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B 稀疏 度 的 算法 和 经 典 SVM 方法 。 文献 [10] 利 用 深度 学 习 特 征 
© 高 光谱 图 像 Chyperspectral image, HIS) 为 每 个 像素 记录 MAARRE ER ELI A2, EH ÉRA P4 A kE 
C. 数 百 个 光谱 带 ， 提 供 高 光谱 分 辩 率 和 增强 的 稳健 图 像 分 类 潜 ”的 图 像 数据 中 提取 深度 特征 ， 使 用 稀疏 表示 分 类 框架 特征 分 
E JUA, ZA, MEER T RRN REJ, a 25. SALLE AD i EI 2/1 p 13 a RI A RES 23 H F 
MÉR RAKAM KA, HE m EEE ”高 光谱 图 像 分 类 ， 能 够 从 未 标记 数据 中 自 适 应 地 学 习 合 适 的 
算法 将 HSI 数据 投影 到 较 低 维 空间 , 包括 传统 主 成 分 分 析 申 ， ”特征 表示 。 
Fisher 线性 判别 分 析 LDADI 以 及 这 些 方法 的 许多 变 体 。 在 文献 [12] 中 ， 已 经 证 明 在 PCA 投影 子 空 间 中 提取 的 二 
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在 文献 [6] 提 出 一 种 基于 高 斯 混合 模型 和 马尔 可 夫 随 机 HÈ Gabor 特征 对 于 HSI 分 类 是 有 效 的 。 在 文献 [13] 中 ， 提 出 
场 的 高 光谱 图 像 分 类 方法 ， 将 局 部 保持 非 负 和 矩阵 分 解 以 及 局 ”了 两 种 基于 Gabor 特征 和 导数 特征 的 融合 分 类 算法 ， 以 实现 
部 Fisher 判别 分 析 部 署 为 预 处 理 ， 以 减少 高 斯 混合 模型 分 类 HSI 分 类 的 优越 性 能 。 文 献 [14] 提 出 一 种 基于 三 维 Gabor 特 
器 的 数据 维 数 ， 同 时 保留 数据 内 的 多 模 态 结构 。 文 献 [7] 提 出 ”” 征 的 协同 表示 方法 用 于 HSI 分 类 ， 该 方法 能 够 降低 分 类 方法 
种 基于 核 的 基于 RBF 核 的 SVM 高 光谱 图 像 分 类 特征 选择 ” 的 计算 复杂 度 。 
方法 ， 可 以 根据 计算 系数 的 大 小 来 实现 特征 小 子 集 和 特征 排 本 文 研 究 了 现 有 高 光谱 图 像 分 类 方法 ， 考 虑 到 原始 光谱 
序 这 两 个 属性 。 由 于 SVM 对 HSI 维 数 不 敏感 ， 而 Markov ”特征 和 空间 特征 ， 提 出 一 种 基于 Gabor 特征 和 决策 融合 的 高 
随机 场 (Markov random field, MRF) 可 以 捕获 复杂 的 空间 背景 ”光谱 图 像 分 类 算法 ， 该 方法 可 以 将 相关 光谱 进行 智能 分 组 ， 
HE, SVM-MRF 已 经 用 于 HIS 的 分 类 。 在 PCA 投影 子 空间 中 提取 每 组 中 的 Gabor 特征 , 并 对 空间 特 
另外 还 有 其 他 关于 HIS 分 类 的 方法 ， 如 局 部 二 值 模式 、 征 进行 降 维 ， 最 后 使 用 决策 融合 规则 对 高 斯 混合 模型 分 类 器 
稀 玻 表示 和 机 器 学 习 。 文 献 [8] 中 采用 局 部 二 值 模式 来 提取 局 ”的 分 类 结果 进行 合并 ， 最 终 实现 高 光谱 图 像 的 分 类 。 
部 图 像 特 征 ， 并 此 采用 特征 级 融合 和 决策 级 融合 局 部 二 值 模 
式 特征 ， 最 后 采用 极限 学 习 机 作为 分 类 器 。 文 献 [9] 提 出 一 种 


T 


DS 


收 稿 日 期 : 2018-09-06; 修 回 日 期 : 2018-10-19 基金 项 目 : 重庆 市 教委 课题 (KJ1729408) ; 重庆 市 教委 教改 重点 项 目 (162071) 
作者 简介 : 杨 秀 杰 (1979-) ， 女 ， 四 川 广安 人 ， 副 教授 ， 硕 士 ， 主 要 研究 方向 为 数字 媒体 应 用 技术 (cqyangxj@sina.com); 高 丽 〈1977-) ， 女 ， 重 庆 
人 ， 讲 师 ， 硕 士 ， 主 要 研究 方向 为 数据 处 理 技术 . 


子 空间 中 的 Gabor 纹理 用 于 提取 空间 光谱 特征 。 
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间 的 维度 。 
GMM) 分 类 器 以 获得 本 地 标签 。 最 
(logarithmic opinion pool, LOGP) 决 
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本 文 分 类 方法 


相关 系数 矩阵 用 于 对 光谱 带 进行 分 组 ， 然 后 PCA 投影 
在 特征 提取 
后 ， 采 用 局 部 保持 非 负 矩阵 分 解 (locality preserving 
-negative matrix factorization, LPNMF) 用 于 减少 特征 子 空 
随后 应 斯 混合 模型 (Gaussian mixture model, 

后 应 用 对 数 意 见 池 

策 融 合 方法 合并 所 有 分 类 
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结果 。 分 类 融合 系统 被 设计 用 于 高 光谱 图 像 分 类 任务 ， 具 体 
流程 如 图 1 所 示 。 
原始 图 像 
Y 
频带 分 组 
组 1 组 2 组 m 
Y Y Y 
GT GT GT 
Y Y Y 
LPNMF LPNMF LPNMF 
Y Y Y 
GMM | GMM GMM 
y 
> 决策 融合 < 
分 类 标签 
图 1 本 文 分 类 方法 流程 图 
Fig. 1 Flow chart of classification method 
1.1 波段 选择 
高 光谱 图 像 被 相关 系数 矩阵 (correlation coefcient matrix, 


CCM) 划 分 为 连续 空间 。 
而 8 Gc») ez HU TR 
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其 中 :ww 是 f(x%,y) 的 平均 值 ， 


Bit Sy) 表示 原始 的 高 光谱 图 像 ， 
真实 图 像 ， 互 相关 函数 定义 为 


STI th yr) us Gey) -nu,] 


SE [2Xtieen-«T | 
«XE eG») 的 平均 值 ，1 表 示 行 


h(Lk)- 


(1) 


大 表示 列 数 。 


1 M N 
u= HELI y) (2) 


(3) 


Band Number 


Band Number 
图 2 相关 系数 矩阵 


Fig.2 Correlation coefficient matrix 
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而 子 空间 之 间 的 统计 依 
(波段 数 为 220) 
留 在 CCM 对 角 线 上 
这 些 带 彼此 高 度 相 
吸水 带 波段 ， 则 去 


分 区 子 空间 提供 增加 的 类 分 离 ， 
us 例如 ， Indian Pines Zt 
自然 被 划分 为 五 个 组 。 这 五 组 对 应 于 驻 
的 五 个 浅 色 块 ， 如 图 2 所 示 。 在 每 组 中 ， 
关 ， 由 于 光谱 带 {104-108} 和 {150-163} 是 
除 这 些 波段 。 
1.2 Gabor 变换 

Gabor 滤波 器 是 由 高 斯 包 络 调 
为 取向 依赖 带 通 滤波 器 ， 并 且 可 
提取 的 物理 结构 的 方向 和 比例 。 然 而 ，HSI 在 很 宽 的 频谱 范 
围 内 包含 大 量 的 频谱 信息 ， 并 且 波 长 间隔 相对 较 小 ， 这 增加 
了 分 类 任务 中 的 处 理 复杂 性 和 统计 病态 调节 。 为 了 提高 效率 ， 
本 文 在 Gabor 特征 提取 之 前 考虑 了 降 维 , 由 于 在 PCA 投影 特 
征 空间 中 提取 的 Gabor 特征 已 被 证 明 可 以 提供 非常 好 的 结 
果 ， 本 文 在 特征 提取 过 程 中 采用 了 这 种 降 维 策略 。 
在 二 维 的 坐标 (ab) H, Gabor 滤波 器 , 包括 实 部 和 虚 部 ， 
J] 以 表示 为 


Ef 


制 的 正弦 函数 ， 可 以 被 视 
以 有 效 地 捕获 HSI 空间 特征 


z| 


a? + y? 
20 2 


ee) 7 


其 中 : 6 代表 正弦 因子 的 波长 ， 儿 是 相位 偏 移 ，7 是 空间 纵横 
比 〈 默 认 值 在 文献 [10] 中 设置 为 0.5) 指定 Gabor 函数 支持 的 
椭圆 度 。w=0 和 w=zx/2 返回 Gabor 滤波 器 的 实 部 和 虚 部 。 


a =acosO+DsinO 


(4.53.0. o. 7) =exp - 


(5) 
b «—asin--bcosO (6) 
在 式 (5)(6) 中 ，6 表 示 Gabor 核 之 间 的 取向 分 离 角 ， 参 数 
的 标准 推导 包 络 线 由 5 和 空间 频率 带宽 bw 确定 。 
ô llIn2 22 +1 
s ie -= E 
在 本 文 提出 的 融合 分 类 系统 中 ， 首 先 考虑 了 每 个 带 通 的 
八 个 方向 c, E, 3m, m 5m 370, ， 并 对 系统 参数 进行 了 优 


848284 8 
化 , 得 到 了 最 优 的 分 类 结果 。, 在 图 3 中 给 出 了 不 同 5 ，pw 和 
PC ( 主 成 分 的 数量 ) 的 总 分 类 精度 (overall classification 


accuracy, OA). 


E 
O XE IBI 


OA (%) 


à z 3 
6 o uw a re uc a) A 3 

S e a a e V ^ 
: w x a PC 


33 不 同 参数 的 整体 分 类 准确 度 

Fig.3 Overall classification accuracy of different parameters 

可 以 观察 到 5 和 PC 的 最 佳 值 分 别 为 65 和 3， 在 所 本 文 
分 类 算法 中 ， 选 择 bw=4 ， 而 该 参数 对 分 类 性 能 不 敏感 。 
1.3 局 部 保持 非 负 和 矩 阵 分 离 

LPNMF 结合 了 非 负 矩阵 分 解 (non-negative matrix 
factorization, NMF) 和 局 部 保持 投影 (locality-preserving 
projection，LPP) 的 优点 ， 这 导致 仅 使 用 加 法 运算 的 基于 部 件 
的 表示 ， 其 中 内 部 几何 结构 也 保留 在 嵌入 空间 中 。LPNMF 
的 优点 可 以 使 其 成 功 对 HIS 进行 降 维 。 另 外 LPNMF 结合 
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GMM 分 类 器 能 够 实现 HIS 的 高 精度 分 类 , LPNME 的 目标 函 “， 据 集 包括 16 个 陆地 覆盖 类 别 ， 其 中 一 些 包含 少量 样本 。 由 于 

数 为 这 样 的 少量 训练 样本 所 获得 的 概率 分 布 不 能 很 好 地 表示 它们 

p 的 统计 特性 ， 所 以 本 文选 择 了 像素 较 多 的 八 类 子 集 ， 对 于 这 

o- Y slog +y |+an & ^ p dus T 

EE Xj 数据 集 ， 平 均 每 组 50 个 训练 样本 和 8624 个 测试 样本 〔 测 

. 试 和 训练 的 比率 约 为 21.6:1) 被 设置 用 于 实验 。 本 文 仿真 实 
参数 ， 式 (8) 的 第 一 部 4 

参数 ， 冻 (8) 的 第 “部 分 是 验 是 在 笔记 本 电脑 的 MATLAB 2013a 中 实现 的 , 笔记 本 电 及 
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其 中 :Y={9j}=UV? ,而 4 是 正 由 


非 负 和 抢 阵 分 解 的 通常 目标 函数 ， 它 使 用 X 和 7 了 之 间 的 配置 为 250G 内 存 和 i7 处 理 器 ，win10 操作 系统 。 
Kullback-Leibler 散 度 。 而 RR 用 于 在 降 维 子 空间 V(t1<4d) 中 的 为 了 量化 所 本 文 算法 的 有 效 性 ， 将 本 文 算法 与 其 他 算法 
点 之 间 强 制 一 个 几何 局 部 约束 。 进行 比较 ， 其 他 方法 包括 : 传统 的 SVM、SVM- 马 尔 可 夫 随 
机 场 模 型 (SVM-Markov Random Field, SYM-MREF)、 局 部 化 的 
z»»0 slog bog% W,, (9) M ; ipd 
js gn Vja Fisher 判别 分 析 - 高 斯 混合 模型 (Local Fisher Discriminant 


Erg. w 是 测量 原始 空间 XxX 中 "m 距离 的 边缘 权重 和 矩阵， Analysis-Gaussian Mixed Mode, LFDA- GMM)、 局 部 保持 非 负 


该 矩阵 遵循 LPP 理论 ， 用 于 保持 数据 分 布 的 固有 几何 形状 。 矩阵 分 解 -高 斯 混合 模型 (locality-preserving non-negative 
下 面 的 乘法 法 则 用 于 最 小 化 函数 上 并 估计 符 阵 U fU V. matrix factorization-GMM, LPNMF-GMM), J5 ili Fisher 判别 
Y Gur) Yo avia) 分 析 - 高 斯 混合 模型 -马尔 可 夫 随 机 场 模型 (LFDA-G-MRF)、 局 
Wig llig EX dO ”部 保持 非 负 和 矩阵 分 解 -高 斯 混合 模型 -马尔 可 夫 随 机 场 模型 
(LPNMF-G-MRF)、 决 策 融 合 高 斯 混合 模型 -局 部 化 的 Fisher 
Tr xa, 1] 判别 分 析 (DG-LFDA) 和 决策 融合 高 斯 混合 模型 -局 部 保持 非 
p» D uma) 4 XB EEAHR(DG-LPMN). Fd 4 给 出 不 同 算法 对 高 光谱 数据 的 
eis x XiaMi a 分 类 结 : R, 
v, se| Zut +A] ie 2 nan) 1) 
Xinlti.g 
day (a Vna) 
其 中 : 是 V 的 tw，7 是 nxn 单 位 矩阵 ，L 是 W 的 图 拉 普 拉 
斯 算 子 。 


1.4 高 斯 混合 模型 和 决策 融合 
GMM 可 以 被 视 为 两 个 或 更 多 个 普通 高 斯 分 布 的 组 合 ， 


在 典型 的 GMM 表示 中 ，X={x%h 在 W 中 的 概率 密度 函数 
被 写 为 上 个 高 斯 分 量 的 总 和 。 


PX- Ý aN Gps) (12) 


(a) 地 面 图 像 (b) SVM 


Ep: N(X. t2) 表示 混合 物 的 必 高 斯 分 量 。Q ,4,3 是 混 
合 物 组 分 的 数量 ， 由 参数 矢量 9={ow, 上 ,2} 表示 。 一 旦 确定 


了 每 GMM 的 最 佳 组 件数 K ， 利 用 迭代 优化 策略 对 混合 模型 
参数 进行 估计 。 
本 文 的 研究 是 在 Gabor 特征 和 诀 策 融合 系统 的 背景 下 进 (c) SVM-MRF (d) LEDA-GMM 


行 的 。 一 个 多 分 类 器 系统 ， idet peius 解决 了 小 
样本 问题 ， 保 证 了 成 组 带 的 鲁 棒 识 别 。 研 究 了 对 数 分 类 ; 
(logarithmic opinion pool, LOGP) 的 分 类 ant. LOGP 需要 使 

每 个 分 类 器 的 个 体 后 验 概率 p(w|x) ，(i=1,2,…,n) 来 估计 全 
类 成 员 关系 : 


| 


3-DIn(wlx 
i=l 


—logC(w;|x) -TI« log p; (w; |x) 
izl 


(e) LPNMF-GMM (f) LFDA-G-MRF 


其 中 : abi 4RASARRUEDSEDAES WEKRE, i ENK 
器 索引 ，7 是 分 类 器 的 数量 ， j 在 分 类 器 组 中 检测 1 类。 
2 ”实验 结果 与 分 析 

本 文 实验 对 象 是 Indian Pines 高 光谱 遥感 数据 集 , 该 数据 
集 是 由 美国 国家 航空 航天 局 红外 成 像 光谱 仪 (administration's 
airborne visible/infrared imaging spectrometer, AVIRIS) 传 感 器 
收集 。 该 数据 集 代表 了 一 个 植被 分 类 场景 ， 具 有 145*145 像 
素 和 220 个 光谱 带 ， 空 间 分 辩 率 为 20m。 原 始 的 印度 松树 数 g) LPNMF-G-MRF (h) DG-LFDA 
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(i) DG-LPNMF 


(j) LFDA-GMM 


图 4 Indian Pines 数据 集 地 面 实 况 和 


Fig.4 Ground truth and classification map of Indian Pines dataset 


从 图 中 可 以 看 


,本 文 算 法 对 高 光谱 医 
其 他 8 种 分 类 算法 ， 这 是 因为 本 文 算法 能 够 对 相 邻 特征 


优 于 
的 局 部 信息 和 高 光谱 区 
效果 明显 优 于 其 他 几 种 分 类 算法 。 


像 分 类 效果 最 好 ， 


像 的 多 模型 结构 进行 保护 ， 故 其 视觉 


表 1 给 出 了 高 光谱 图 像 的 总 体 分 类 准确 度 (overall 
classification accuracy, OA) 和 Kappa(k) 系 数 的 分 类 结果 , 本 次 


实验 取 20 次 实验 了 


表 1 


均值 。 
不 同 算法 的 高 光谱 


像 分 类 性 能 指标 


Table 1 Performance indices of hyperspectral image classification 


using different algorithms 


算法 OA(%) Kappa 
SVM 77.46(+1.71) 0.7337(+0.0195) 
SVM-MRF 88.04(*1.14) 0.8578(50.0135) 
LFDA-GMM 64.17(€2.04) 0.5767 (50.0218) 
LPNMF-GMM 73.01(+1.67) 0.6806(+0.0194) 
LFDA-G-MRF 82.09(+2.89) 0.7840(+0.0351) 
LPNMF-G-MRF 88.17(+2.59) 0.8491 (+0.0303) 
DG-LFDA 93.23(+1.15) 0.9192(+0.0170) 
DG-LPNMF 95.60(*0.74) 0.9473 (+0.0089) 
本 文 算 法 99.22(+0.46) 0.9907(+0.0055) 


AR 1 中 数据 可 以 得 至 


比 LPNMF-GMM 高 0.3101 


生 能 优 于 其 他 8 


OA, Kappa 和 标准 差 都 是 最 优 怕 


PF 现 有 算法 。 图 5 给 


4， 本 文 算法 的 OA 达到 99.22%， 

比 传统 SVM 高 约 22%, Eb SVM-MRF 高 11% , 比 DG-LPNMF 
高 4%; 本 文 算法 Kappa 达到 0.9907, EE LEDA-GMM 高 0.414， 
， 比 SVM 高 0.257， 本 文 方法 的 
能 ， 对 于 高 光谱 图 像 的 分 类 


准确 率 的 关系 


OA 在 不 同 训练 样 
清晰 、 更 平滑 的 分 类 医 


3 ”结束 语 


本 文 提出 一 种 高 
[高 斯 混合 模型 的 融合 


EY 


从 图 5 可 以 看 
即 随 着 训练 样本 的 增加 ， 本 文 算法 各 
本 数量 对 LFDA-G-MRF 算法 的 影 


首 加 ，LFDA-G-MRF 的 OA 值 增幅 较 


的 数量 下 都 优 了 
， 说 明 本 文 算法 的 


6 谱 图 像 分 类 算法 ， 
策 算 法 ， 解 决 传 


I 


E 的 问题 。 


6 谱 特 征 而 忽略 空间 


矩阵 对 HSI 的 相关 频 乔 


H RE 


和 了 训练 样本 量 和 分 类 


,训练 样本 数量 对 本 文 算 法 的 影响 较 小 ， 
J OA 增幅 较 小 ， 训 练 样 
向 最 大 ， 随 着 训练 样本 的 
大 。 另 外 ， 本 文 算法 的 
他 算法 ， 从 而 导致 更 


有 效 性 。 


即 基 于 Gaber 特征 


统 分 类 算法 只 考虑 


tF 


该 算法 


首先 使 用 相关 系数 


行 分 组 ， 然 后 使 


Gabor 纹理 从 每 


杨 秀 杰 ， 等 : 基于 降 维 Gabor 特征 和 决策 融合 的 高 光谱 图 像 分 类 


日 中 提取 空间 频谱 特征 。 为 了 保留 相 邻 像素 的 局 部 结构 ， 
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N 


入 维特 征 执行 GMM 分 类 ， 调 用 基于 LOGP 的 决策 融合 方法 
将 这 些 结果 组 合成 每 个 像素 的 单个 分 类 决策 ,实验 结果 表明 ， 
本 文 算法 优 于 传统 算法 和 现 有 先进 算法 ， 且 分 类 准确 率 在 很 
少 样本 数 下 就 能 达到 较 高 的 OA 值 ， 能 够 产生 更 清晰 、 平 滑 


的 分 类 图 。 


OA(%) 


LPNMF 将 高 维 子 空 间 投 影 到 较 低 维 子 空间 上 ， 之 后 对 


GI 


- 本 


95 


90 


一 6 一 本 文 算法 
—IL-DG-LPNMF 


-——DG-LFDA 


—m-LPNMF-G-MRF 


= 一世 一 LFDA-G-MRF 


40 45 50 55 60 
训练 样本 数量 
图 5 分 类 准确 率 与 训练 样本 量 的 关系 


Fig.5 Relationship between classification accuracy and training 


sample size 
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